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Small-scale locomotion plays an important role in biology. Different modelling ap- 
proaches have been proposed in the past. The simplest model is an infinite inexten- 
sible two-dimensional waving sheet, originally introduced by Taylor, which serves 
as an idealized geometrical model for both spermatozoa locomotion and ciliary 
transport in Stokes flow. Here we complement classic steady-state calculations by 
deriving the transient low-Reynolds number swimming speed of such a waving sheet 
when starting from rest (small-amplitude initial- value problem) . We also determine 
the transient fluid flow in the 'pumping' setup where the sheet is not free to move 
but instead generates a net fluid flow around it. The time scales for these two prob- 
lems, which in general govern transient effects in transport and locomotion at low 
Reynolds numbers, are also derived using physical arguments. 

Keywords: swimming problem; low Reynolds number; locomotion; 
micro-organisms; cell motility 

1. Introduction 

The locomotion of microorganisms plays a vital role in biology. Examples include 
the locomotion of mammalian spermatozoa during reproduction, or the swimming 
of bacteria and algae to locate better nutrient sources (Childress 1981; Bray 2000). 
Microorganisms adopt different swimming strategies from those exploited by larger 
animals. For humans, fish, or birds, swimming and flying are usually accomplished 
by imparting momentum to the fluid. This strategy, however, no longer works in 
the world of microorganisms, where inertia plays little role and viscous damping is 
paramount. The Reynolds number, a dimensionlcss parameter characterizing the 
ratio of inertial to viscous forces in the surrounding fluid, ranges from 10~ 6 for 
flagellated bacteria to 10 -2 for spermatozoa (Brennen & Winet 1977). Many of such 
small organisms propel themselves by propagating progressive waves along their 
flagclla. The geometrical characteristics of these microorganisms and the waves they 
propagate have been reviewed by Brennen & Winet (1977). For example, for animal 
spermatozoa, typical wavelengths range from 11 — 65 (mi and wave amplitudes range 
from 4—15 /im, while the ratio between their swimming speeds and the wave speeds 
range from 0.07 to 0.3. Taylor (1951) initiated the studies on hydrodynamics of 
microorganisms by modelling flagclla swimming. After his pioneering work, much 
progress has been made towards addressing the basic mechanisms of propulsion for 
swimming microorganisms, and we refer to Brennen & Winet (1977) and Childress 
(1981) for detailed reviews. In this paper, we focus on unsteady effects in low- 
Reynolds number locomotion. We start below with a brief overview of relevant 
works before summarizing the approach and outline of our paper. 
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In 1951, Taylor first showed that self-propulsion of a waving sheet is possible in 
a viscous fluid in the absence of inertia. The sheet moves in a direction opposite to 
that of the wave propagation, and the steady state swimming speed was computed 
explicitly to be a 2 k 2 c/2 at leading order in the wave amplitude, where a, k and c are 
the wave amplitude, wave number and wave propagation speed respectively. Drum- 
mond (1966) extended Taylor's results to larger amplitudes of motion. Reynolds 
(1965) and Tuck (1968) included inertial effects in the analysis and found that the 
swimming speed decreases with Reynolds number. In the inertial realm, Wu (1961, 
1971) also considered the problem of a waving sheet, but with finite chord and in the 
inviscid limit, as a model for fish propulsion. Recently, Childress (2008) discussed 
the nature of the high Reynolds number limit of Taylor's swimming problem and 
applied the results to the idea of recoil swimming, where propulsion is achieved by 
movements of the center of mass and center of volume of the body. 

Despite its simplicity, the waving sheet serves as an idealized geometrical model 
for both spermatozoa locomotion and ciliary transport. Cilia are short flagella beat- 
ing collaboratively to produce fluid motion. They are important in many biological 
transport processes, for example, the transport of mucus in the respiratory tract of 
humans. Blake (1971) represented the envelope formed by the cilia tips as an infi- 
nite impenetrable oscillating surface, assuming that the cilia are sufficiently closely 
packed together. The use of an infinite sheet to model finite length organisms is 
particularly well justified for elongated and flat organisms such as Paramecium or 
Opalina (Blake 1971). If the waving sheet is not allowed to move, there will be 
a net flow of the fluid in one direction. In this case, the sheet acts as a pump for 
transporting the fluid. We shall refer this as the 'pumping problem' in this paper. In 
addition, when a second waving sheet is present above the original waving surface, 
one recovers the problem of peristaltic pumping, a useful fluid-transport mechanism 
in physiology and many industrial processes. Readers are referred to reviews on this 
subject (Jaffrin & Shapiro 1971; Fauci & Dillon 2006). 

Another scenario of interest is the propulsion of microorganisms near solid 
boundaries. Natural microorganisms often swim in narrow passages such as the 
swimming of spermatozoa in the cervix (Suarez & Pacey 2006). In addition, during 
most laboratory examinations, the presence of coverslips imposes solid boundaries 
near the microorganisms. Reynolds (1965) adopted Taylor's swimming sheet model 
to study swimming near solid walls for small-amplitude waving motion; another ap- 
proximation, the long wavelength limit, was considered by Shack & Lardner (1974) 
and Katz (1974). The oscillating wall models by Smclscr et al. (1974) and Shukla 
et al. (1988) and the layered fluid medium model by Shukla et al. (1978) were also 
proposed to study the interaction between the cervical wall and the cell. 

The fact that many biological fluids are non-Newtonian has also received consid- 
erable attention. Chaudhury (1979) first extended Taylor's swimming problem to 
viscoelastic fluids. The same problem was then considered by Sturges (1981) using 
a more rigorous integral constitutive equation. Fulford et al. (1998) modified the 
resistive force theory to model the swimming of a spermatozoon in a general lin- 
ear viscoelastic fluid. Recently, Lauga (2007) revisited Taylor's original calculation 
using more realistic, nonlinear non-Newtonian fluid models. 

This brief literature review shows that most studies in small-scale biological 
locomotion focus on solving for the swimming speed of a model organism. Since 
previous work derived solutions for steady-state swimming only, we propose in this 
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paper to go beyond the steady limit and study a prototypical time- varying situation, 
namely the initial-value problem of a model microorganism (Taylor's waving sheet) 
starting from rest. Physically, such a process is governed by the small time scale 
necessary for vorticity created at the swimmer surface to propagate diffusively into 
the fluid, and belongs to the general class of unsteady Stokes problems. 

As expected, this transient swimming process is also dependent upon the de- 
velopment of the propagating wave from rest; we present here a general analytical 
treatment for a waving sheet which develops its frequency (or phase speed) from 
rest to the steady state in an arbitrary manner. An analytical formula describing 
the general transient propulsion speed of the sheet is then derived, complement- 
ing Taylor's well-known steady state solution. We also solve the transient pumping 
problem, where the sheet is not free to move, but instead entrains the surrounding 
fluid in an unsteady manner. Unlike their steady counterparts, these two problems 
are not equivalent in the unsteady case because of time-dependent inertial forces. 

The paper is structured as follows. In §2, the swimming sheet problem is mathe- 
matically formulated with the appropriate non-dimensionalization, governing equa- 
tions, and boundary conditions. In §3, the calculations at first and second order 
are presented. The results of §3 are then applied to study the transient pump- 
ing problem (§4) and the swimming problem (§5). Finally, a physical discussion of 
our results and a derivation of the time scales involved in transient low-Reynolds 
number swimming is offered in §6. 

2. Formulation 

We consider here an infinite sheet swimming in an incompressible fluid, similar 
to the model proposed by Taylor (1951). We also allow the wave of displacement 
along the sheet to include not only normal but also tangential motion (Blake 1971; 
Childress 1981). Here, the waving motion is observed in the frame moving at the 
unknown swimming speed. As the swimming speed is time-dependent during the 
transient motion, the reference frame is non-incrtial and hence a fictitious (inertial) 
force has to be introduced (see below). In the moving frame, the position of material 
points, {x m ,Hm)i on the waving sheet is written as 

x m (x, t) = x + a A cos(kx — io(i)t — </>), 
y m {x,t)=bAsm(kx-w(t)t), (2.1) 

where a and b are dimcnsionless, A is a typical wave amplitude and 4> is the phase 
difference between the longitudinal and transverse motion. The angular frequency 
u>(t) is an arbitrary function of time which initially starts from zero and eventually 
reaches a steady state value of Woo- The functions describe a traveling wave with 
tangential motion of wavelength A = 27r/fc moving in the positive x direction at 
a time-varying speed c = ui(t)/k. We keep oj(t) arbitrary in the analysis below to 
derive general formulas describing the transient net flow and transient swimming 
velocity. Specific examples for ui(i) will then be considered for illustration. The 
transverse traveling wave considered by Taylor (1951) is obtained when a = 0, and 
the b = case represents a longitudinal traveling wave. The unknown swimming 
speed of the sheet is denoted by U{t) and is assumed to occur in the direction 
opposite to that of the wave propagation (sec figure [T] for notation) . 
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Figure 1. Two-dimensional swimming sheet and notation. 

(a) Non-dimensionalization 

We non-dimcnsionalizc lengths by l/k and velocities by the wave speed c. The 
time scale in this problem shall characterize changes in velocity in the fluid dur- 
ing the transient motion. The appropriate, physically-motivated time scale, is that 
required for the vorticity created by the start-up of the sheet to propagate dif- 
fusively into the fluid over the characteristic length scale, and therefore time is 
non-dimcnsionalizcd by l/k 2 v, where v is the kinematic viscosity of the fluid. For 
a wavelength on the order of 10 /xm, and in water, the time scale considered is on 
the order of milliseconds. The angular frequency, w(t), is non-dimcnsionalized by 
its steady state value, uJoq. The dimcnsionless quantities (starred) are summarized 
as follows, with p being the density of the fluid, 

x* = — t* = 1 u* = - p* = P uj* = — (2 2) 

l/k l/(fc 2 i/)' c' pvck 1 Woo' 

and the Reynolds number is given by Re = c/vk = Lo OQ /vk' 1 . The dimensionless 
position of material points on the sheet is written as 

x m *(x*,t*) = x* +eacos(x* - Re uj*(t*)t* - 0), 

y m *(x*,t*) = ebsm{x* - Re u*(t*)t*), (2.3) 

where e = Ah. We assume that the wave amplitude is much smaller than the 
wavelength, and derive the results in the limit where K<1. 



(b) Governing Equation 

We consider an incompressible Newtonian fluid surrounding the sheet. Since we 
have a two-dimensional setup, the continuity equation, V • u = 0, is satisfied by 
introducing the stream function ij}(x,y,t) such that u = dtp/dy and v = —dip/dx, 
where u = ue x + ve y . In the laboratory frame, the governing equation is the Navier- 
Stokes equation 

/ 9^+ /9 (u-V)u = -Vp + ^V 2 u, (2.4) 

where p is the pressure and fi the dynamic viscosity of the fluid. In the non-incrtial 
frame considered in this paper, as we accelerate with the swimming sheet, a time 
dependent uniform fictitious force F(t) = (F x (t), 0, 0) has to be introduced into the 
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equation as 

p— +p(u- V)u = -Vp + ^V 2 u + F. (2.5) 

Since the swimming sheet accelerates uniformly in the x-direction, the non-zero 
component of the fictitious force, F x , is uniform (independent of space) and is a 
function of time only. The non-dimcnsionalizcd form of the equation (with non- 
dimcnsionalization described above) is given by 

du* 

— +Re (u* • V*) u* = -V V + V* 2 u* + F*. (2.6) 

dt* 

For small-scale biological locomotion, we consider the low Reynolds number limit 
Re — > 0, where the convective term vanishes, resulting in the unsteady Stokes' 
equation 

^- = -V>* + V* 2 u* + F*. (2.7) 

dt* v ' 

Upon taking the curl of the equation, both the pressure gradient and the fictitious 
force terms vanish, resulting in an equation for the ^-component of the vorticity, 
lu*. With the relation u* = — S7* 2 ip* , the equation for the stream function is given 
by 

d 



V* 2 V* 2 ^* = 0. (2.8) 
dt* J 

Hereafter, we shall mostly deal with the dimensionless quantities, and therefore the 
stars will be omitted for simplicity. 

(c) Boundary Conditions 

The unknown swimming velocity of the sheet is denoted by —U(t)e x . In the 
frame moving with the swimming sheet, the velocity of the fluid in the far field 
(y = oo) is therefore given by U(t)e x . Hence, the far-field boundary conditions are 

|(x m ,oo)= U(t)e y . (2.9) 

On the swimming sheet, the boundary conditions are given by the velocity compo- 
nents of a particle of the sheet 

\(x y )= e f{t) [bcos (x — Re ■ ui(t)t) e x + asin (x — Re ■ w(t)t — (f>) e y ] , 

(2.10) 

where f(t) is a function of time, defined by f(t) — d(u(t)t) /dt. The conditions 
simplify to 

V ^ \(x m ,y m )= ef(~t)[b cos (x)e x + a sin (x - <f>)e y ] , (2.11) 

in the low Reynolds number limit (Re — ► 0). Note that Re = LOr^jvk 2 can be 
re- written as Re = [l/i//c 2 ]/[l/<x' 00 ], and can therefore be interpreted as the ratio 
between the relevant time scale for viscous diffusion into the fluid and the typical 
time scale of the wave (its period). The Re — > limit physically means that the 
diffusion of vorticity occurs much faster than the propagation of the wave, and 
therefore the wave appears to be stationary on the time scale where viscous diffusion 
is taking place. 
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3. Analysis 

In this section, we seek regular perturbation expansions for the stream function, 
swimming speed and the pressure in powers of e in the form 

{V, U,p} = e{^, 17l.pi} + e 2 {^2, U 2 ,p 2 } + .... (3.1) 

The order e and order e 2 solutions are presented in the following subsections. 

(a) First-order solution 
At order e, the governing equation is given by 

^-V 2 )vVl = 0. (3.2) 

Expanding the boundary conditions on the sheet using Taylor expansion, they be- 
come at order e 

i>i,y\(x,o) = a f{t)sm(x -0), (3.3) 
V>i,x l(x,o) = bf(t) cosx, (3.4) 

lpl,y l(x,oo) = U X (t), (3.5) 

ipi,x \(x,oc) = 0. (3.6) 

The initial condition is that the vorticity in the fluid (u> = — V 2 ^) is initially zero, 
i.e. uji | t=0 = (-V 2 tpi) | t=0 = 0. 

To allow an easy implementation of the initial condition, we solve the prob- 
lem using the Laplace transform method. A similar technique has been applied to 
study the transient solution of Stokes' second problem (Erdogan 2000). The Laplace 
transform of the stream function ipi is defined by the relation 

/>oo 

4> 1 (x,y,s)= Mx,y,t)e- St dt, (3.7) 
Jo 

where s is the Laplace variable and tilde variables represent transformed quantities. 
Taking the Laplace transform of the governing equation at this order, equation (|3 . 2[) 
becomes 

(s - V 2 ) V 2 ^i = 0. (3.8) 

The boundary conditions in Laplace domain are given by 

i>i,y lo,o) = af(s)sm(x~(p), (3.9) 

V>i,x |(x,o) = bf(s)cosx, (3.10) 

i>X,V l(x,oo) = U^s), (3.11) 

i>l, x |(x,oo) - 0. (3.12) 

The solution satisfying all the boundary condition is found to be 

_ / e ~Vsy\ f(s) > 

■4'i = Ui\y+ — ■=- + = ' asin0(e-^+ Ty - e" y )cosa; 

\ V s \/s + 1 - 1 



V7+T- l 



(b^/sTT + acos(f>)e- y - (b + acos^e-^^v sin a;. (3.13) 
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(b) Second-order solution 

We proceed to the analysis at order e 2 using a similar procedure. The governing 
equation at this order is 

I- - V 2 ) VV 2 = 0, (3.14) 



dt 



with boundary conditions 



/I I A\ d 2 ^ 1 I J, • I IO 1 

V2, y \(x,o) = -acos{x - <P)-q^ \y=o - bsmx -g^T \y=o, (3-15) 
W2.x \{x,o) = -acos{x - 9)-g-2- \y=o ~ bsmx g x g y ^ =0 ' (3.16) 

^2,y \(x,oo) = U 2 (t), (3.17) 
~4>2,x \(x,oo) = 0, (3.18) 

where we have used Taylor expansion to obtain the boundary conditions on the 

moving sheet surface. Using Laplace transform, the governing equation for the 
stream function becomes 

(s - V 2 ) V 2 -0 2 = 0. (3.19) 
The boundary conditions are now 



i>2, y |(x,0) = /(s) 



b 2 \/s + 1 sin 2 x + ab — ■ since sinfx — 6) — a 2 cos 2 (a; 

y/a + T-1 



(3.20) 

4>2,x |( X ,0) = 0, (3.21) 
i>2,y \{x,oo) = U 2 ( S ), (3.22) 
$2,x \(x,oo) = 0. (3.23) 
The solution satisfying all the boundary conditions at this order is obtained to be 

i>2 = U 2 \y+ ) - (b 2 J7Tl -a 2 + ab^JL cos A e"^ 



q 2 sin2(ft + afe S sinc/> ) (e^^ 7 ^ - c~ 2y ) sin 2x 

Vs + 1 - 1 / 



2(Vi + 4-2) 

- r Zi___ ( a 2 cos 20 + 6 2 V7TT + a6^=L ) ( c -^+^ - e~ 2y ) cos 2x. 

2( % /JT4-2)V VF+I-l/ 



(3.24) 



4. Pumping Problem 



The analysis of §3 can first be applied to the pumping problem, where the sheet 
is not allowed to move and its oscillatory motion entrains a net fluid flow along 
the cc-direction. Since the sheet does not move, we are in the laboratory frame of 
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reference and the introduction of the fictitious force is unnecessary in this case. 
The velocity in the far field is zero for any finite time, which is analogous to that 
in Stokes' problems; therefore, U\ = U2 = 0, and we then take the x-average of the 
horizontal velocity (denoted by (...)), leading to 

(u 1 )(y,s) = ^ = 0, (4.1) 
dy 

(u 2 )(y,s) ^^M^M (V-v/sTl -a 2 + ab—=^= — - cos <j\ e"^. (4.2) 

As expected from the e — > — e symmetry, no net flow occurs at order e. At order e 2 , 
there is a net fluid flow and the result of equation (|4.2|) describes the (dimcnsionlcss) 
transient velocity of the driven flow in the Laplace domain. For any finite distance y, 
as time goes to infinity, the average horizontal velocity, by the final value theorem, 
asymptotes to 

(«2>(y,*-> 00) = i(fo 2 + 2a6cos0-a 2 ), (4.3) 

which is the appropriate steady-state value (Blake 1971; Childress 1981). The di- 
rection of the net fluid flow is governed by the oscillation mode of the swimming 
sheet. For a transverse wave (a = 0), there is a net fluid flow in the positive x 
direction, which is the direction of the wave propagation. For a longitudinal wave 
(b = 0), the fluid flows in a direction opposite to the wave propagation. If there 
is a combination between transverse and longitudinal motion (a ^ 0,6 / 0), the 
direction of fluid flow is governed by the competition between the amplitudes of 
the transverse and longitudinal motion and their phase difference, as described by 
equation (|4.3[) . Notably, there are values of a, b and cf> that will yield a zero net fluid 
flow, where the propulsion caused by the longitudinal motion exactly balances that 
produced by the transverse motion. 

The detail evolution of the transient flow depends on the development of the 
propagating wave, described by the function cu(t) (or f(t) = d (u>(t)t) /dt). Examples 
will now be given to illustrate the use of equation (|4.2p to determine the transient 
velocity of particular driven flows in the time domain. 



(a) Example 1: tu(t) = 2 arctan (t/T) /n 

Here, we consider the angular frequency of the wave which develops according 
to the function oj(t) = 2 arctan(i/T)/7r, where T characterizes the time taken for 
the angular frequency to reach its steady state value, non-dimensionalized by the 
viscous diffusion time scale. When T is equal to zero, it represents the case where 
the wave starts impulsively. Computing the function f(t) = d(u>(t)t) / dt and its 
Laplace transform f(s), the transient velocity of the driven flow in Laplace domain 
follows from equation (|4.2p as 

, , . , 2Ci(sT) [sin(sT) - sTcos(sT)] + [cos(sT) + sTsin(sT)] [tt - 2Si(sT)l 

«2 ){y,s) = 

Z1TS 

b 2 Vs+T- a 2 +ab . L cos^)e-^ y , (4.4) 

v s + 1 — 1 / 
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(«2) («2) (U2) 




10= 10' 10° 10' 10 2 10» 10" 10 a 10' 10" 10' 10 2 10 3 10* 10= 10' 10" 10' 10 2 1 3 1 0* 

t t t 

(a) Transverse wave (b) Longitudinal wave (c) Combined wave 

Figure 2. Average dimensionless horizontal velocity in the fluid for the pumping problem, 
(m2)/moo, at the non-dimensional position y = 1, as a function of the dimensionless time, 
in the case where the wave frequency increases according to u(t) = 2 arctan(f/T)/7r. 
(Semi- log plot) (a): transverse wave (a = 0); (b): longitudinal wave (6 = 0); (c): combined 
transverse and longitudinal wave (a = 6 and <j> = 7r/4). In each graph, from left to right: 
T = 0.01, 0.1, 1, 10, 100. 



where Ci and Si are the cosine and sine integrals, defined respectively by 

«(*)=- r^A, si(*)= r^-dt. (4.5) 

The Laplace transform inversion can be easily implemented numerically (Valko 
& Abate 2004) . The evolution of the transient swimming speed in the time domain is 
computed at the non-dimensional position y = 1 and shown in figure [2] for different 
cases: transverse wave (figure [2^), longitudinal wave (figure Eb) , and the combined 
wave where a = 6 and <fi = 7r/4 (figure Ek) . The parameter T is varied from 0.01 to 
100 in each case, and as expected an increase in T leads to an increase in the time 
necessary for the flow speed to reach its steady state value. 



(6) Example 2: Impulsive motion 

We consider now the extreme case where the angular frequency of the propagat- 
ing wave attains its steady state value instantaneously, i. e. if T is the dimensionless 
start-up time, we have T — 0. The solution to that problem exists as a regular 
perturbation only in the case where 6 = 0, i.e. when the wave has no transverse 
amplitude. This is due to the fact that in the limit t — * 0, and as in Stokes' first 
problem for the impulsive motion of a plate in a viscous fluid, an infinite shear is 
initially created along the surface of the body. Since a Taylor expansion in the wave 
amplitude is used to obtain the second-order boundary conditions, equation (|3 . 1 5|) . 
the presence of an infinite shear leads to an infinite boundary condition unless 6 = 0, 
in which case the Taylor expansion does no evaluation into the fluid domain (see 
§6 for further discussion). However, similarly to Stokes' first problem, the solution 
is well-behaved when 6 = 0, which is the case we now consider. In this case, the 
function f(t) is a unit step function and its Laplace transform is f(s) = 1/s. Again, 
as expected from the e — > — e symmetry, no net fluid flow occurs at order e. Net flow 
occurs at order e 2 . By equation (|4.2p . the x-average horizontal velocity in Laplace 
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(u 2 ) (u 2 ) 




t ' y 

(a) (b) 

Figure 3. (a): Average dimensionless horizontal velocity in the fluid for the pumping prob- 
lem, (u2)/uoo, at the non-dimensional position y = 1, as a function of the dimensionless 
time for an impulsive longitudinal wave 6 = 0. (Inset: same plot on a semi-log scale); (b): 
Evolution of the average dimensionless horizontal velocity in the fluid in the y direction 
at different values of dimensionless time for an impulsive longitudinal wave, from left to 
right: 4 = 0.1,1,5,10,15. 



domain in this case is given by 



which we inverse Laplace transform, yielding 

(« a >(tf,t) = -yerfc^) I (4.7) 

where erfc is the complementary error function. The fluid flow occurs in the direc- 
tion opposite to the wave propagation. The result of equation (|4.7|) describes the 
(dimensionless) transient velocity of the driven flow, and is illustrated in figure [3J 
The vorticity perturbation created by the instantaneous motion of the sheet prop- 
agates diffusively into the fluid, y ~ y 7 !, with a final result very similar to that of 
Stokes' first problem. For any finite distance y, as time goes to infinity, the average 
horizontal velocity asymptotes to — a 2 /2. In other words, the entire fluid will even- 
tually be pumped to move with an average horizontal velocity of — a 2 /2 for very 
large times. In dimensional form, the average horizontal velocity is given by 

, . . , a 2 A 2 k 2 c „ / y \ . , , 

<«>(**) = — erfc (^fc) ( 4 - 8 ) 

at leading order in sheet amplitude. 



5. Swimming Problem 

Next, we consider the swimming problem, where the sheet is free to move. We 
are now in a frame moving with the sheet, a non-inertial frame of reference (§2 b). 
In this problem, the swimming speed of the sheet is yet to be determined. An 
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additional condition is required to determine U\(s) and ^(s) in equations (|3.13[) 
and (|3.24p respectively. For this purpose, Newton's second law will be applied on 
the sheet in the x-direction. By the periodicity of the problem, we consider the 
forces (per unit width) acting on one wavelength of the sheet. The forces are: (i) 
the horizontal force per unit width acted on the sheet by the fluid (denoted Ffl u id); 
and (ii) the fictitious force per unit width due to the accelerating reference frame 
(denoted Factitious)- In the frame moving with the sheet, its acceleration is zero, 
and therefore Newton's second law on the sheet under this non-inertial frame has 
the form Ffl u id + -Factitious = 0. The magnitude of the fictitious force, -Factitious, is 
given by the mass per unit width of the sheet times the acceleration of the reference 
frame. The fictitious force is acting in a direction opposite to the acceleration of the 
moving frame, so that -Factitious = (mass of the sheet per unit width) dU / dt. Below, 
we expand both Factitious and F<fl u id in powers of the small parameter e, and enforce 
Newton's second law at each order. 



(a) First-order solution 
At order e, and in dimensional variables, F^uid is given by 

r-2-n p2ir 

pvc I (ui y + v\ x ) |y=o dx = pvc I C~ 1 (U\\fs)dx = 2-KpvcCT 1 (jj\\fs) 1 (5.1) 
Jo Jo 

where C~ l is the inverse Laplace transform operator. The fictitious force, Factitious, 
at this order reads (p s h\ck 2 i/)dUi/dt, where p s and h are the density and thickness 
of the swimming sheet respectively. Applying Newton's second law, we have 

2itpvcC- 1 (U 1 *f~s) + p s hXck 2 u^- = 0. (5.2) 

dt 

Taking the Laplace transform of the equation yields 

2irpvcU 1 yfs + p a h\ck 2 vsUi = 0. (5.3) 

Solving for Ui(s), we have Ui(s) = 0, which implies Ui(t) = 0. As expected, no 
self-propulsion occurs at order e. We therefore proceed to the analysis at order e 2 . 



(6) Second-order solution 

A similar analysis is undertaken at this order, with the difference that the order 
e pressure, pi, will be required for the calculation of the fluid force, F^uid, at order 
e 2 . The pressure is found by integrating over x the Navier-Stokes equation at order 
e in the horizontal direction. The forces, Factitious and Ffluid, are expanded on the 
sheet using Taylor expansion. At order e 2 , Ffl u ;d is given by 

pvc I [—(2ui, x — pi)bcosx + b8mx{ui,y + vi, x ) y +U2 tV +V2,x]\y=tidx 
Jo 



2irpvcC 



U 2 \fs- -^-^3 (b 2 \/s + 1 - a 2 + ab cos cj)(\/s + 1 - -y/s + 1)) 



(5.4) 
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The fictitious force, -Factitious, at this order reads (p s h\ck 2 v)dU2/dt. Again, applying 
Newton's second law on the sheet leads to 



' (&Vs + 1 - a 2 + ab cos 4>(Vs + 1 - \fs + 1)) 



+pji\ck 2 v^- = 0. (5.5) 



Taking the Laplace transform, we have 



Un 



y/s (b 2 y/s + 1 - a 2 + abcos </>(y/s + 1 - y/s + 1)) 



/As 



p (1/fc) 



st/ 2 = 0. 



(5.6) 



We define the dimensionlcss parameter M. = (p s / p)(hk), which is the product of 
the density ratio of the sheet to the fluid and the ratio of the thickness of the sheet 
to the characteristic length scale (wavelength). The swimming velocity at order e 2 
is then determined, in the Laplace domain, as 

/» 



^(*) = 'V 7T (b 2 V7TT - a 2 + abcos^V^TT - x/^s + I)) ■ (5.7) 
Z{MyJs + 1) 

The steady state swimming velocity, by the final value theorem, is given by 



1 



Uoo = - (b 2 + 2ab cos <j) - 



(5.8) 



which agrees with previous steady-state results (Blake 1971; Childress 1981). When 
a = 0, it corresponds to Taylor's result of a transverse traveling wave. The swimming 
sheet propels itself in a direction opposite to the wave propagation. The 6 = limit 
corresponds to the case of a longitudinal traveling wave, and the swimming sheet 
propels in the direction of wave propagation. If there is a combination of transverse 
and longitudinal motion, similar to the pumping problem, the swimming direction 
is governed by the competition of the amplitudes of the transverse and longitudinal 
motion and their phase difference, as described by equation (|5.8p . Similarly to the 
pumping problem, there are values of a, b and <p that will produce a zero propulsion 
speed. The transient swimming velocity of the sheet is given by equation (|5.7|) in 
the Laplace domain for arbitrary function f(s). Here again, we consider different 
examples of /(s) and compute the transient swimming velocity in the time domain. 



(c) Example 1: ui(t) = 2arctan(i/T) /n 

Here, we consider the angular frequency of the wave which develops according 
to the function uj(t) = 2 arctan (t/T) /it, the same function as for the pumping 
problem above. Therefore, the functions f(t) and /(s) are the same as in §4a. By 
equation (|5.7p . hence, the swimming velocity in Laplace domain is given by 

- 2Ci(sT) [sin(sT) - sTcos(sT)] + [cos(sT) + sTsin(sT)] [tt - 2Si(sT)] 

2(S) ~ 2ns (AV5TT) 

x (bVs + 1 - a 2 + ab cos 0(Vs +T - y/s + 1)) . (5.9) 
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u 2 u 2 u 2 




t t t 

(a) Transverse wave (b) Longitudinal wave (c) Combined wave 

Figure 4. Dimensionless swimming speed of the sheet, U2/U00, as a function of the di- 
mensionless time, in the case where the dimensionless wave frequency evolves according to 
ui(t) = 2 arctan(t/T)/7r, and for A4 — 1 and different values of the wave start-up time T. 
(Semi-log plot), (a): transverse wave (a = 0); (b): longitudinal wave (b — 0); (c): combined 
transverse and longitudinal wave (a = b and cf> = tt/4). In each graph, from left to right: 
T = 0.01, 0.1, 1, 10, 100. 

The transient swimming velocity in the time domain is obtained by performing a 
Laplace transform inversion numerically. The results are displayed in figures [Hand [5] 
for different modes of oscillation: transverse wave (figures HJt and [5ji) , longitudinal 
wave (figures HJd andO)), and the combined transverse and longitudinal wave where 
a = b and <fi = 7r/4 (figures and [St). In figure [H the mass ratio is fixed to be 
A4 = 1 and the wave start-up time T is varied between 0.01 to 100. In figure [5j 
we fix T = 1, while the parameter M. is varied between 0.01 to 100; the limit 
where M = is also computed and shown as dashed lines in figure CD This limit 
corresponds to the case where the swimming sheet is mass-less, which is the relevant 
limit for biological organisms. 



u 2 u 2 u 2 




10"* 10 s 10" W 10" 10 s 10" 10 1 " 10- 4 10- 2 10° 10 2 10" 10« 10" 10'° 10-' 10-* 10° 10 s 10" 10 6 10» 10" 

t t t 



(a) Transverse wave (b) Longitudinal wave (c) Combined wave 

Figure 5. Dimensionless swimming speed of the sheet, U2/U00, as a function of the di- 
mensionless time, in the case where the dimensionless wave frequency evolves according to 
ui(t) — 2 arctan(f/T)/7r, and for T = 1 and different values of the mass ratio M. (Semi-log 
plot), (a): transverse wave (a = 0); (b): longitudinal wave (6 = 0); (c): combined trans- 
verse and longitudinal wave (a = b and 4> = 7t/4). In each graph, from left to right: 
M = 0.01, 0.1, 1, 10, 100. The dashed lines (the left-most in each plot) refer to the case of 
a mass- less sheet i.e. M = 0. 
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10 2 10° 10 z 10" 

t 

(a) ui(t) = 2arctan(t/T)/7r 



t 

(b) u{t) = l-e- t/T 



Figure 6. Dimensionless swimming speed of the sheet, U2/U00, as a function of the dimen- 
sionless time, for T = 1, J\A = 1, and for two different different wave start-up function 
w(t). (Semi-log plot). Blue (solid line): transverse wave (a = 0); Red (dashed line): lon- 
gitudinal wave (6 = 0); Green (dotted line): combined transverse and longitudinal wave 
(a — b and <f> = 7r/4). 

(d) Example 2: w(t) = 1 - exp(-i/T) 

To illustrate the difference between different types of wave start-up, we now 
consider the case where the angular frequency increases exponentially as: uj(t) — 
1 — exp(— t/T). Again, computing the corresponding function f(t) and its Laplace 
transform f(s), the swimming velocity in Laplace domain, by equation (|5.7[) . reads 



(1 + 2sT) 



2a(l + My/s)(l + sTf 



l-a A 



ab cos(/>(vi 



1 - 



!))• 
(5.10) 



Here again, we invert the Laplace transforms numerically. The influence of the di- 
mensionless parameters T and M. on the transient swimming velocity show trends 
similar to the first example above. The detailed evolution of the swimming velocity 
is however, and as expected, quantitatively different due to the different time devel- 
opment of the angular frequency. The two different evolution profiles arc illustrated 
in figure O for the case where T = 1 and M. = 1, and for the three different waves. 

In the limit where the swimmer has no mass, M. — 0, the Laplace transform 
inversion is simplified and analytical formulas for swimming velocity are obtained 
for all the three different modes of oscillation. The swimming velocity for a mass-less 
sheet, for a transverse wave, reads 



U 2 (t) = b 2 



f V< erf (Vt) 



(2t~3T~2tT + 2T 2 ) erf 



t(T-l) 



2T^ 



47^3/2^7^- 



whercas for a longitudinal wave, it is 

2 



(5.11) 



(5.12) 
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for combined motion where, for example, a = b and <j) — 7r/4, the formula is 



U 2 (t) = 2 -±^a> 



V2 



e-tVt erf(Vt) «"* (* - 3T - 2T + 2I«) erf ( 
+ 2 



4T 3 / 2 VT T7 T 



4 

V2 - 2 



* (2t 3 / 2 - 4TVt + T 3 / 2 erfi 

2T 3 7v! 



3-* I --1 



(5.13) 



where erf is the error function and erfi is the imaginary error function. 

(e) Example 3: Impulsive motion 

Similarly to the pumping problem, the case of impulsive motion is singular 
except in the case where 6 = for which the impulsive swimming speed is well 
behaved. In this case, and as in the pumping problem, the function f(t) is a unit 
step function and its Laplace transformation is f(s) — 1/s. Hence, the swimming 
velocity of the sheet, in Laplace domain, is given by 



& 00 



2s {M^ + 1) 



(5.14) 



The inverse Laplace transform of equation (|5.14j) yields the leading order swimming 
velocity in the time domain 



U 2 (t) 



a 

Y 



1 w) erfc 



(5.15) 



As in the corresponding pumping problem, the transient swimming velocity involves 
the complementary error function, erfc(x), as illustrated in figurc[7J For large times, 
the velocity asymptotes to — a 2 /2. In dimensional form, the swimming velocity of 
the sheet is given by 



1 — exp 



VP 

h 2 P 2 s 



t erfc 



hp s 



A 2 k 2 c , 



(5.16) 



at the leading order in the wave amplitude. 



6. Discussion 

In this paper, we have studied two different unsteady Stokes flow problems, and 
obtained explicit analytical formulas in the Laplace domain for their transient mo- 
tion. In the pumping problem, the fluid is driven by the waving motion of a fixed 
sheet, and the Laplace transform of the leading-order dimensionlcss average hori- 
zontal velocity is described by equation (|4.2| . In the swimming problem, where the 
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Figure 7. Dimensionless swimming speed, U2/U00, for the impulsive longitudinal wave, as 
a function of the reduced non-dimensionalized time, t/M 2 (Inset: same plot on a semi-log 
scale) . 

sheet is free to move, the leading-order swimming velocity of the sheet is given by 
equation (|5.7|) . Higher order terms may be obtained in a similar fashion. 

We first note the difference between the form of equation ()4.2|) and equa- 
tion (|5.7p . In the steady limit, the pumping and swimming problems differ only 
by a change of reference frame, so the fixed waving sheet pumps a uniform amount 
of fluid at a velocity equal to minus the steady swimming speed of the free-swimming 
sheet. In the unsteady case however, the sheet is accelerating, and therefore is sub- 
ject to additional forces. As a result, the final formula for the swimming speed 
involves the sheet mass (through the dimensionless parameter A4). 

To obtain the solution in the time domain, Laplace transform inversion of equa- 
tion (|4.2p and equation (|5.7p is required. Analytical formulas have been obtained 
for several cases where the Laplace transform inversion is simple. It can be further 
noted that, if a mass-less waving sheet (M = 0) is considered in the swimming 
problem, the Laplace transform inversion is particularly straightforward for the 
case of longitudinal wave (6 = 0). By equation (|5.7|) . we inverse Laplace transform 
to yield the simple formula U(t) = —e 2 a 2 f(t)/2. 

For more complex cases (M. 7^ 0, 6 7^ 0), the Laplace transform inversion is 
easily implemented numerically for any admissible /(£) or f(s). Particular examples 
have been given for illustration. We have introduced the dimensionless time T 
characterizing the time scale over which the angular frequency reaches its steady 
state. As illustrated in figure [5] for the pumping problem, the development of the 
flow field can be divided into two stages. In the initial stage, the development of 
the transient velocity depends on the transient motion of the waving sheet, and 
the relevant time scale is T: the smaller the values of T, the shorter the time is 
required to reach a given velocity. For large times, the waving motion of the sheet 
has effectively reached its steady state, and the subsequent development of the flow 
field is dominated by the viscous diffusion alone. The values of the parameter T 
are no longer relevant, and the curves for different values of T collapse onto the 
same envelop for large times (see figure [2]). In addition, and as could be expected, 
the detailed evolution of the net velocity depends on the different waving modes 
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of the sheet. Similar trends are observed in the swimming problem, as displayed 
in figure |U In that case, since the sheet is free to move, its mass — characterized 
by the dimensionless parameter M. — comes into play in its transient swimming 
velocity, and leads to the presence of a third relevant time scale (see below). As 
illustrated in figure [SJ the mass of the sheet dictates its acceleration in the initial 
stage, and the sheet with the smallest mass accelerates the fastest. 

The case where T = 0, i.e. where the wave starts impulsively, has also been 
discussed. We find that, within the framework of the perturbation expansion pro- 
posed in this paper, this singular start-up behavior for the sheet leads to a singular 
pumping and swimming solution — except for the case where the wave motion is 
purely longitudinal, i.e. b = 0. Physically, as the wave is impulsively starting, it cre- 
ates an initially infinite shear rate immediately above the sheet. Since the boundary 
condition for the second order flow is found by the "sampling" of the first order 
flow by the first-order shape of the sheet (Taylor expansion), any case where the 
sheet protrudes into the fluid (b ^ 0) leads to singular boundary conditions for 
the second-order flow. This singular behavior could be resolved by incorporating 
the advective inertial terms in the Navier-Stokes equation. For example, in the case 
where u>(t) = 1 — exp(— t/T), it can be shown that, for a wave with a non-zero trans- 
verse amplitude, the convective term scales as ~ 1/T at the sheet position (y = 0). 
Hence, for the convective term to be negligible, it is required that Re/T <C 1 or 
Re <C T, which means that for the results of the present work to be uniformly 
valid for all time, a transverse wave needs a finite start-up time. For a longitudinal 
wave however, the convective term remains order unity for all time, even when the 
wave is propagated instantaneously, and therefore the impulsive motion problem is 
well-posed. 

A close examination of the results obtained in the impulsive case for the longitu- 
dinal wave allows us to understand physically the origin of the time scales involved 
in the unsteady pumping and swimming processes. In the pumping problem, the net 
velocity of the fluid averaged over one wavelength occurs in the direction opposite 
to the wave propagation, with magnitude given by equation (|4.8|) . The scaling for 
the time evolution of the fluid velocity in that case is straightforward and similar 
to that of Stokes' first problem. It is a diffusive scaling, and the fluid below a dif- 
fusive front propagating as y ~ \f~\A into the fluid is pumped roughly at the steady 
velocity. In the swimming problem, the sheet moves in the same direction as the 
longitudinal wave with a transient propulsion speed given by equation (|5.16|) . In 
that case, the scaling for the relevant time scale for start-up of the swimmer arises 
from a consideration of the balance between the inertial force of the sheet and the 
shear stresses exerted by the fluid. The inertial force is given by the mass per unit 
width and per unit wavelength of the sheet times its acceleration. Over one wave- 
length, the small-amplitude swimmer has a mass per unit width equal to p s h, and 
the acceleration scales as U s /t, where U s denotes the sheet velocity. On the other 
hand, the shear stress on the sheet exerted by the fluid is on the order of pUf/y, 
where Uf denotes the fluid velocity and y ~ \f~\A is the typical size in the y direction 
of velocity gradients in the fluid. Balancing shear stresses in the fluid with inertia 
in the sheet leads to p s hU s /t ~ /j,Uf/y/id. At steady state, we have U s ~ Uf, and 
the balance suggests a typical time scale, £m, given by tM ~ P 2 h 2 /p 2 v, which is the 
time scale involved in equation (|5. 16|) . Compared to the typical diffusive time scale, 
tfj ~ l/k 2 v, we have tuj^D = M 2 - For typical swimming microorganisms, the 
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Figure 8. Evolution of the dimensionless swimming speed, U2/U00, for the case of longitu- 
dinal wave and M = 1, given by the result of this paper (solid lines) and the quasi-steady 
approximation (dashed lines) for different values of T as indicated in the figure. 



ratio of swimmer density to fluid density is about one, and the ratio of thickness h 
to wavelength 1/k is on the order of 0.01 (Brennen & Winet 1977; Childress 1981), 
so the steady state swimming problem is reached much earlier than the steady-state 
pumping problem. 

Finally, the results of this paper can be compared with those given by a quasi- 
steady approximation, where the steady-state formulas of Taylor are extended to the 
unsteady case by replacing the value of the steady-state frequency by the instan- 
taneous value oj(t). For illustration, we consider the profile cu(t) = 2 arctan (t/T) /tt 
in the case of longitudinal wave, with M. = 1. The comparison is displayed in fig- 
ure [SI the solid lines show the evolution of the propulsion speed obtained by our 
analysis while the dashed lines are obtained using the quasi-steady approximation. 
Two values of T are employed for illustration. For T = 100, the development of 
the wave frequency is slow compared with viscous diffusion, and the quasi-steady 
approximation leads to a reasonable agreement with our results; the inertial effects 
are relatively unimportant in this limit. In the other limit for T = 0.01, the wave fre- 
quency increases quickly compared to the viscous time scale, and the quasi-steady 
analysis significantly over-estimates the swimming speed; the inclusion of inertial 
effects as carried out in this paper is thus critical in this case. 

In conclusion, our results complement Taylor's classical swimming sheet calcu- 
lation by deriving the transient pumping and swimming motion of the sheet. Gen- 
erally, the two time scales derived above, which control the startup of the flow sur- 
rounding the sheet (tp ~ \jk 2 v), and the startup of the swimmer (iju ~ M 2 /k 2 v), 
are expected to govern transient effects in transport and locomotion at low Reynolds 
numbers. In addition, our study could be extended to the case of viscoelastic fluids, 
for which the relaxation time scale of fluid becomes important, and more complex 
transient processes such as the switching of rotation direction in bacterial flagclla. 

This work was funded in part by the US National Science Foundation (grants CTS-0624830 
and CBET-0746285 to Eric Lauga). 
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